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Abstract: We extend Choptuik's scaling phenomenon found in general relativistic critical 
gravitational collapse of a massless scalar field to higher dimensions. We find that in the 
range 4 < D < 11 the behavior is qualitatively similar to that discovered by Choptuik. In 
each dimension we obtain numerically the universal numbers associated with the critical 
collapse: the scaling exponent 7 and the echoing period A. The behavior of these numbers 
with increasing dimension seems to indicate that 7 reaches a maximum and A a minimum 
value around 11 < D < 13. These results and their relation to the black hole-black string 
system are discussed. 
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1. Introduction 

General Relativity (GR) in higher dimensions (D > 4) has been receiving increasing atten- 
tion in recent years. This is largely motivated by the phenomenological "brane- world" and 
"large-extra-dimensions" scenarios that are rooted in String theory, which is intrinsically 
higher-dimensional, involving 10 or 11 dimensions. But regardless, GR in itself does not 
favor the 4D case especially; rather it is independent of the dimension, and so without any 
additional reasons the dimension of space-time may and should be considered a parameter 
of the theory. 

Adopting this point of view, one finds that certain properties and solutions of GR 
change qualitatively in higher dimensions. In particular, black hole uniqueness breaks 
down in D > 4. This implies the existence of black objects with non-spherical horizon 
topologies. Examples include the caged black hole - black string system in a space-time 
with compact extra dimensions, see 0] for a review, and the rotating black hole - black 
ring configurations ||, ||]. The phase transitions between different horizon topologies seem 
to lead to a compromise of cosmic censorship in higher dimensions and to be accompanied 
by energy outbursts (see e.g. Q.) Moreover, there exist critical dimensions above which 
the qualitative properties of the spoken phase transition change H ](| [7| . 

In this paper we study the critical collapse of a spherically symmetric massless scalar 
field in D-dimensional spacetime, in search of similar non-trivial dimensional effects. The 
critical behavior of black hole formation is a striking example of the surprising phenomena 
to be found in gravitational physics. It is understood that in the dynamical gravitational 
collapse of matter fields, initial configurations with very low density will usually disperse 
to infinity, while extremely dense clumps will collapse to form a black hole. What happens 

1 Another example of the qualitative dependence of GR on dimension is the chaotic behavior of the 
spacetime near a spacelike singularity, discovered by Belinskii, Khalatnikov & Lifshitz (BKL). Specifically, 
above the critical dimension Dbkl = 10, the system becomes non-chaotic, see the recent review fej. 
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in the intermediate case between these two extremities is much less obvious. In 1992, 
Choptuik H studied this limiting case numerically, modeling the collapsing matter as a 
spherically symmetric configuration of massless scalar field coupled to gravity. For each 
family of initial data (e.g. gaussian, lorentzian etc.), parameterized by an amplitude p, there 
is a critical amplitude p*. If p > p* the final state is a black hole (supercritical collapse), 
otherwise the field disperses leaving behind empty flat space (subcritical collapse). The 
following discoveries were made regarding the behavior of this system near the threshold 
amplitude p* in AD: 



In supercritical collapse arbitrarily small black holes are formed as p — > p*. In 
this limit the black hole mass M scales as Moc (p — p*) 7 , with the pre- factor and p* 
characteristic of the specific family of initial data. The exponent, however, is universal 
and independent of the shape of initial data. Its numerical value is 7 ~ 0.374. 
Later it was also discovered that in subcritical collapse the maximal scalar curvature 
encountered before the field disperses scales as R ma x on (p* — p)~ 27 [10|. 



The critical solution itself is universal. Namely, following a short transient stage, the 
spacetime converges to a universal solution independent of the initial data, remains 
such for a while and then either collapses or disperses, depending on whether p > p* 
or not. 



• The critical solution is discretely self-similar. Designating the critical solution (col- 
lectively for the scalar field and the metric) by Z*(r,t), this means that Z#(r,t) = 
Z*(r e^,te^), with the empirically found A ~ 3.44. The field and metric functions 
pulsate periodically with ever decreasing temporal and spatial scales, until a singu- 
larity is formed at the accumulation point, r = t = 0. 



Similar results were also observed for other kinds of matter and extensive literature 
has been written on the subject. For a review see 11] and the references therein 2 . 



There are several reasons why we find it interesting to explore the critical phenomena 
in higher dimensions: (i) The mysterious universal numbers that appear in critical collapse, 
such as the scaling exponent and the echoing period, are known only empirically. One may 
wonder if it is possible to use the dimension as a probe that can provide insights or hints 
into the origin of these numbers; (ii) The Choptuik critical phenomena, or more precisely 
its time-symmetric version, is related to the phase transition in the black-string-black- 



hole system [14]. It might be possible to learn about the behavior of that system near the 
merger point where the topology changes from that of a black string to that of a black 
hole by better understanding the higher-dimensional critical collapse scenario; (iii) The 
massless scalar field collapse is probably the simplest model of dynamical collapse in GR. 
It is worthwhile to understand what kind of new features one can expect to find there in 
higher dimensions. 



2 While a majority of authors do not go beyond 4D, Garfinkle et. al. |12| have obtained 7 and A in 6D 
and Birukou et. al. p^], got 7 in the critical scalar field collapse in 5D and 6D. 
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To realize these goals we designed a numerical code that evolves the space-time and 
scalar field equations and allows us to study the critical solution in different dimensions. As 
is often the case in numerical analysis, we encountered unexpected problems that demanded 
special treatment. In addition to the familiar methods used in 4D like "constrained" 
evolution and mesh refinement, the code uses analytical series expansion near the origin 
and "synthetic friction" (or "smoothing") to handle the severe divergences appearing in 
higher dimensions. We employed this code for 4 < D < 11. 

In all our simulations we find that the near-critical collapse proceeds in the discretely 
self-similar (DSS) manner as in the AD case. For each dimension we obtain the (logarithmic) 
period of the pulsations, A, and the scaling exponent, 7, which we define such that the 
dimension of length is [length] oc \p — p*| 7 . These results are summarized in table [l]. 
In the examined range of dimensions we find that A decreases and 7 increases 3 with D. 
While we could not check this explicitly, since our numerical tools did not take us further 
than D = 11, we were tempted to extrapolate our results beyond that. This extrapolation 
seems to indicate that A reaches a minimum and 7 has a maximum somewhere in between 
11 < D < 13, see figures |l(] and |ll|. In any case, the dimensional dependence of these 
variables is smooth (apparently lacking in divergences). Moreover, 7 varies slowly with 
respect to D, and the critical exponent for mass in D dimensions is 7 maS s = (D — 3)7, so as 
the dimension grows the mass of a black hole forming above the threshold increases steeply 
with p — . 

In the next section we formulate our problem: we derive the equations and define the 
variables. Then in section |3| we describe our numerical scheme. The results are presented 
in section [4]. These and their relation to the black-hole-black-string system is discussed in 
section [|. In the same section we state some open questions. 

2. Equations and variables 

We turn now to formulate the physical problem at hand. In spherical symmetry, we describe 
the D dimensional asymptotically flat space-time in double-null coordinates (u,v) by the 
metric: 

ds 2 = —a(u, v) 2 dudv + r(u, v) 2 dti 2 D _ 2 , (2-1) 

where dfl 2 j_ 2 is the metric on a unit S^* -2 sphere and the axis r = is chosen to be where 
u = v. 

The action of the massless scalar field (j) minimally coupled to gravity in D dimensions 
is given by 

S = — / R^d D x - [ g ab d a ^d b (p^d D x. (2.2) 



16vrG u 

The Einstein equations, R a b = SuGcf)^^ derived from this action under the metric 
24]) are 

r >uv + (D - 3) '— = 0, (2.3) 



3 There is some decrease in 7 for D = 11 but we can not definitely determine if this reflects a physical 
behavior or caused by the escalation of numerical errors in higher dimensions. 
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where the first two are hyperbolic equations of motion, and the other two are constraints. 

Variation of the action with respect to <p yields the scalar field equation of motion, 
□0 = 0, which expands to 

<p,uv H t, — = °- ( 2 - 7 ) 

2 r 

Following Hamade&Stewart |ll| we formulate the problem as a set of coupled first order 
differential equations. Designating s = \/4irG4> we define 

Dl: w = s u D2: z = s v D3 : f = r u DA : g = r v D5 : d = (2.8) 

a 

Then the complete set of evolution equations is 

El: f tV + (D -3)fs±*£ = 0, 

r 



£2: (P-2)( D - 3 )/ 9 + cV4 +ii ,^ 0| 

2 



a „ 2 o 

CI: /, u - 2^/ + — — rV = 0, 
a D — 2 

C2: g }V - 2dg + ^-^rz 2 = 0, 

51 : z u H = 0, 

2 r 

£-2 /z + gW 

52 : u; „ H = 0. 

2 r 

Some useful scalar quantities include the Ricci scalar curvature, which is given by 

R = - 8 -^ (2.9) 
and the proper time of an observer at the axis, 

T{u) = [ a(u',u')du. (2.10) 



The supercritical collapse of a spherical distribution of matter in D asymptotically flat 
dimensions results in the formation of a D-dimensional Schwarzschild-Tangherlini black 
hole, whose metric reads [[ilfj 

ds 2 = -f{r)dt 2 + /(r^dr 2 +r 2 dn 2 D _ 2 , 

f(r) = 1 - (r /r) D -\ (2.11) 
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The ADM mass of the black hole is M = (D - 2)A D - 2 ro~ 3 /(16vrG) in terms of its 
Schwarzschild radius, Tq. Ad-2 = 2tt d / 2 ~ 1 /T(D/2 — 1) is the area of a unit S D ~ 2 sphere 
and G is the Z)-dimensional Newton constant. 

We define the critical exponent, 7, such that \p* — p| 7 has dimensions of length. Then 
it follows [l^] that in D-dimensions the maximal curvature (having the dimension of 
inverse length squared) achieved in a subcritical collapse scales as R ma x (p* —p)~ 21 and 
the mass of the black hole, forming in a supercritical collapse, behaves as (p — p*) 7mass 
where 7 mass = (D — 3)7. 

Initial data problem, gauge and boundary conditions 

We specify the initial scalar field profile along an outgoing null surface u = Ui = const, 
which we choose to be ui = 0. To complete formulation of the initial value problem we note 
that the choice of metric (|2.1| ) is only unique up to the redefinitions v' = x( v )> u' = C( u )- 
In order to fix this residual gauge freedom we choose the area coordinate, r, along the 
initial null surface u = as r = v/2 4 . In addition, we set a = 1 at the axis on m = v = 0. 
From here we can obtain all the other functions on the initial hyper-surface by integration 
from the origin. 

Equations El — E2 and 51 — S2 are singular at the axis. The physical solution is, 
however, perfectly regular there. This enforces the following boundary conditions on the 
scalar field and the metric functions along the axis r = 

9 = - f = 
w = z, 
d r s = 0, 

d r a = 0. (2.12) 

The source terms which are singular at the axis are then evaluated by applying l'Hospital's 
rule. 

The actual shape we choose for the initial configuration of the scalar field is a gaussian 

shell 



s(v, u = 0) = p exp 



. 2 
v — Ve- 



er 



(2.13) 



where v c and a are constants and the amplitude p is the aforementioned strength parameter 
of the initial data. 



3. The numerical scheme 

Having laid out the framework, we proceed to describe the numerical solution of the above 
equations. We construct a grid in the u-v space as described in figure |l[ The primitive 
computational cell is square, with grid spacings h v = h u = h. The initial ingoing wave is 

4 This is chosen to conform with the conventional definition of characteristic coordinates in flat space, 
v = t + r and u = t — r. 
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Figure 1: The domain of integration. At any moment, in addition to the current outgoing hy- 
persurface, Lq we keep two preceding levels L\ and L,2- Boundary conditions involving d/dr are 
implemented using 3-point derivatives along the shown diagonal line. Mesh refinement is illustrated 
in the topmost rows. The smoothing of z (or d) near the axis at a point marked by a cross is done 
using the values of z (or d) at points on the past light cone of that point marked with circles. 



specified on an outgoing hyper-surface of constant u, which is stored as a horizontal line in 
our grid, see figure]]] . Supposing we know the solution along an outgoing hyper-surface u—h 
then d and z are propagated to u using the equations E2 and SI correspondingly. Then 
we integrate the equations 52, El, C2, D2, DA and D5 from the axis outward along u to 
get w, /, g, s, r and a respectively. Note that we use what is called "constrained evolution" , 
namely the constraint C2 is used in the evolution. This algorithm is iterated until either 
the domain of integration is entirely covered (in the subcritical case) or a black hole forms 
(in a supercritical collapse). The remaining equations are not used directly (but are used 
to determine the boundary conditions) in this scheme, but they must be satisfied by the 
solution. We monitor the equations D3 and CI and verify that they are indeed satisfied 
to a prescribed precision during the evolution. 

We use the 4th order Runge-Kutta (RK4) algorithm to solve all equations. At every 
moment in addition to the solution along the line currently being solved {Lq at u) we keep 
in memory the solution on the two preceding levels: L\ at u — h and L2 at u — 2h. Without 
being too memory consuming, this makes the RK4 algorithm more elegantly applicable 
and also helps us to implement the axis boundary conditions as sketched in figure 0. The 
size of the domain is chosen such that the v max = u max point is located just beyond the 
accumulation point. 
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Series expansion and smoothing near the axis of evil 

The basic scheme described above works well for collapse in D = 4, 5, but in higher 
dimensions the code becomes unstable. Moreover, for D > 5 this scheme is unstable even in 
flat space, where only the wave equations 51 — 52 are solved (with constant g = —f = 1/2). 
This instability clearly arises near the axis where r is small and the discretization errors 
in evaluating the sources are amplified. To illustrate what happens near the origin let us 
consider the wave equation 52 in flat space. Say we are solving for w along a u = const 
ray at the grid-point, labeled 1, next to the axis, labeled 0. First we rewrite the wave 
equation as w tV -\- [{D — 2)/2]i/jr jt) /r — \{,D — 2)/4]z/r, where we used the flat space value 
/ = -1/2. Its formal solution is w(vi) = [(D - 2)/4]/r D / 2 ~ 1 Q r D l 2 ~ 2 z(v')dv' '. Using e.g. 
trapezoidal rule we get 

wi - w = — - 2 —zi - z + 0(h 3 ), (3.1) 



where we utilized the boundary condition ( p. 12 ), wq = zq. Alternatively, approximating the 



integrand to linear order in h by Taylor series about we evaluate the integral analytically 
to arrive at 

w 1 -w = — ^{ Zl -zo) + 0{h 2 ), (3.2) 
Comparing ( |3.1| ) and (|3.2j ) we learn that as D increases the first term in the r.h.s of 



(3.1) (where r\ = h/2 ) grows linearly and dominates while the r.h.s. in the analytical 
formula ( |3.2j ) barely varies. Hence the discrete version (3.1) overestimates w\ by a factor 5 



proportional to D, which clearly causes problems in higher dimensions. 

We employ the series expansion method around the axis in equations 52 and -El 6 
such that the sources in these equations are evaluated at the axis. Then the propagation 
of z and d in the u-direction is executed as before but with w and / obtained from the 
series expansion in all evaluations of sources, required by the Runge-Kutta algorithm. We 
determine the number of points near the origin where the series expansion is required 
empirically for every dimension. We found that optimally this number is approximately 
D points. (For example in D = 8 and for the initial grid-spacing of 10 -3 this is done at 9 
grid points.) 

Unfortunately it turns out that while the series expansion prolongs the life-time of 
the code before crashing, it is not enough to render it completely stable. After trying 
several other methods we ended up by adding an effective "friction" or "smoothing" in the 
^-propagated equations (E2 and 51). Specifically, the value of the function z (or d) at a 
point near the axis is mixed with values extrapolated from neighboring points on the past 
light cone of that point (see figure |l|.) Symbolically, the value of some function y at the point 
marked by a cross in figure |l| is updated according to y x = y e + 3to Yli=i Vi /(I + 
where y e is the value obtained from the evolution equation at the crossed point, and 
yi(u) = 2y(u — h) — y(u — 2h) are the extrapolated values along the 3 directions lying in 



J Methods other than the trapezoidal rule can give another factor, which is still proportional to D, so 
for large D's the difference between the discrete and the analytical formulae increases anyway. 
6 The series expansion in El also averts error amplification in E2. 
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the past light cone, designated by circles on the figure, and u are the weights. The weights 
given to the values obtained from the evolution and from the extrapolation depend on D, 
but they are essentially comparable. (For the above D = 8 example the weights ~ 1.) This 
procedure is in principle equivalent to adding a diffusive term in the evolution equations. 
The smoothing is applied at the same mesh-points as the series expansion. 

This series-smoothing symbiosis is a stable method that satisfies the constraints D3 
and CI. Its convergence rate, however, is less than quadratic but never below the linear. 
This is not a serious problem though. The typical run-time of our code on a PC-class 
computer with satisfactory precision is on the order of minutes in spite of this relatively 
slow convergence. We use that series-smoothing method to obtain all the results reported 
in this paper. 

Mesh refinement 

Close to the threshold amplitude the solution becomes DSS, i.e. it is self-replicating 
on decreasing scales. In order to resolve this solution the grid spacing must be smaller than 
the smallest feature in the solution. A static grid with the spacing necessary for resolving 
several echoing periods will require a large amount of memory and computation time. This 
is very inefficient, and moreover unnecessary, because in the early stages the scales involved 
are still relatively large, so a very dense grid is superfluous to resolve them. This makes 
dynamical grid refinement imperative for a realistically feasible scheme. While there exist 
completely general adaptive mesh refinement techniques |i~7|| , we use a far simpler method 
proposed by Garfinkle 

The scheme relies on the fact that when u is incremented the grid point that lies on 
the axis (at (v = u)) enters the region v < u and effectively leaves the causal past of the 
domain (see fig. ||). Physically the ingoing ray is reflected from the origin and becomes an 
outgoing ray, but grid points which were assigned to it are lost. Thus the active part of 
the grid becomes smaller at large u, exactly when more resolution is needed for resolving 
the small features of the critical collapse. This unhappy situation can be prevented by 
interpolating the remaining points at v > u back into the original array of points, restoring 
the original resolution. Following Garfinkle we do this only when half of the grid points 
cross the axis. Since d = hoii the axis, we get a linear increase in the grid density and a 
linear decrease in h with u. 

We use the code to calculate the spacetime near the critical solution in dimensions 
4 < D < 11. For each D we empirically locate the threshold amplitude p* by a simple 
binary search until p* is found with the desired accuracy, which is typically one part in 
10 10 - 10 15 (the lowest is for D = 11). A typical grid has initially between 2000 and 8000 
grid points in the initial outgoing surface. After several mesh refinements the step size h, 
initially of the order of 10~ 3 to 10~ 4 , is reduced by an order of magnitude or two. 

4. Results 

In all dimensions that we examined, the near critical collapse of the massless scalar field 
proceeds in a manner qualitatively very similar to the familiar AD critical collapse. For p < 
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Figure 2: D = 6: Contours of the scalar field profile in slightly subcritical collapse. After a short 
transient, the field oscillates, approaches the accumulation point where the curvature is maximal 
and then disperses. 

p* the curvature along the axis grows, reaching some maximal value and then diminishes. 
In the strong curvature region the scalar field shows echoing on decreasing scales and 
subsequently disperses, see figure In the supercritical collapse the field again rings but 
the curvature finally diverges and a black hole forms (the formation of an apparent horizon 
is signaled by e.g. the change of sign of g, which indicate that outgoing null rays are tilted 
back to smaller radii and do not escape to infinity.) 

In subcritical collapse we can define the "accumulation point" to be the location of 
maximal curvature 7 . We also label the proper time ( |2.10| ) of an on-axis observer at the 
accumulation point by T* . In figure [|| we plot the scalar field on the axis as a function of 
log(T* — T). After an initial transient the field becomes periodic in log(T* — T), which is 
the attribute of DSS behavior in T. We were typically able to observe about 4 full ringings 
of the scalar field in AD and about 3 for higher D. Prom this data the oscillation period is 
computed, and so we get A(D). 

Naturally, not only the scalar field is periodic, but other metric functions exhibit 
oscillations as well. In figure ^| we plot a at the axis. It shows twice the number of 
oscillation of the scalar field, and decreases sharply on approaching the accumulation point. 
The evolution of the scalar curvature ( |2.9D along the axis is shown in figure [| to oscillate 
with the same period as the scalar field. Since all variables are periodic, the quantity A 
can be derived from any one of them. We found that A's obtained from the various metric 

7 For the exactly critical solution it coincides with the singularity. 
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Figure 3: D = 5: The scalar field on the axis as a function of log(T* — T). The period of osculations 
is A ~ 3.19. The actual data is designated by points. The distance between the points increases 
close to T* indicating the decrease of resolution. 
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Figure 4: D = 5: The metric function a on the axis decays fast as the accumulation point T* is 
approached. Its evolution is accompanied by oscillations whose period (~ 1.6) is half the period of 
the scalar field (~ 3.19). 



functions are consistent, see table [T] for the values of A in the verified dimensions. 

We evaluated the critical exponent 7 from the behavior of the maximal curvature R ma x 
in subcritical collapse, and not directly from the black hole mass scaling in supercritical 
collapse. The reason is that it is easier numerically to follow subcritical collapse than to 
estimate the black hole mass in supercritical collapse, because of difficulties in determining 
the position of the apparent horizon. (In addition, as we demonstrate below, it becomes 
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Figure 5: D = 7: Behavior of the curvature as the accumulation point T* is approached. The 
evolution of the curvature, like other metric functions, is accompanied by oscillations. After each 
pulsation log(l — R) increases by A. The period of the last six echoes is approximately constant 
and equal to A/2 ~ 1.41. Similar behavior (with a different echoing period) is observed for other 
Z?'s as well. 



increasingly difficult to find near-critical black holes in higher dimensions.) 

For each dimension we examined the maximal curvature on the axis in several subcrit- 
ical simulations. The result for D = 4 is plotted in figure |6] and the results for other D's 
are summarized in figures [7| and ||. It is evident from these plots that in all dimensions the 
maximal curvature has a dominant power law scaling in p* — p, with the exponent given 
by the slope of the linear fit in the figures. The critical exponent, 7, is then the minus half 
of value of the slope 8 . 7's for different dimensions are listed in table ffl. 

In addition to the dominant power-law scaling there are some residual "wiggles" or 
fine-structure in the curvature's behavior, as can be seen in the figures. This was predicted 
for DSS solutions and explored in Ji"9| , p0[ where it was concluded that the "wiggling" 
period is A/{2^). To illustrate this, we subtract from the measured R m ax the dominant 
power-law dependence and plot the result for D = 8 in figure ||. By reflection about the 
horizontal axis the data is reorganized to have the doubled period. In all dimensions we 
find the "wiggling" period to agree well with the theoretical prediction, A/7. This is a 
non-trivial test for the numerically computed A and 7. Alternatively, since usually we can 
get A with somewhat higher precision, 7 can be computed from the "wiggles" period. This 
is what we do in the D = 11 case because in this case the linear fitting approach does not 
yield a very accurate result, see figure S. 



^Recall, we define 7 such that \p — p, | 7 has dimensions of length. 
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• Data 




log(p-p) 

Figure 6: D = 4: The maximal curvature on the axis as a function of log(p* — p). The slope of 
the linear fit yields 7 = 0.372 which agrees well with the value cited in the literature (~ 0.374). 



D 


A 


7 


4 


3.37 ± 2% 


0.372 ± 1% 


5 


3.19 ± 2% 


0.408 ± 2% 


6 


3.01 ± 2% 


0.422 ± 2% 


7 


2.83 ± 2% 


0.429 ± 2% 


8 


2.70 ± 3% 


0.436 ± 2% 


9 


2.61 ± 3% 


0.442 ± 2% 


10 


2.55 ± 3% 


0.447 ± 3% 


11 


2.51 ± 3% 


0.44 ± 3% 



Table 1: The echoing-period of the scalar field A and the scaling exponent 7 in different dimensions. 
The error in A represents the variation about the mean A measured on several periods of oscillation. 
The error in 7 originates both from the linear fitting and form the actual numerical errors in the 
measured scalar curvature. However it has been verified by convergence analysis that the dominant 
error is in the fitting. The values of A and 7 that we find for D — 4,6 agree well with the numbers 
found in the literature. 

5. Discussion 

Let us now discuss the dimensional dependence of A and 7. We find that in the examined 
range of D's the echoing period A(D) is a decreasing function of dimension, as shown in 
figure [h]. The observed dimensional dependence is well behaved in the sense that nothing in 
this plot forebodes that A will suddenly blow up. Assuming that this regularity continues 
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Figure 7: The logarithm of the maximal curvature on the axis as a function of log(p» —p) for differ- 
ent p's in various dimensions. The slope of the linear fit yields —27. Some "wiggles", characteristic 
of the DSS solutions, occur about the linear fit. 



(otherwise, a discontinuity in A will mark a phase transition) we add, in the same figure, 
a 2 nd order polynomial fitting and continue it slightly beyond the last data point. This 
extrapolation indicates that A might have a minimum 9 at about D ~ 12. 

Our computed j(D) increases up to D < 10 but for D = 11 it appears to decrease, 
although the error-bar is consistent with a constant or even slightly increasing value, (see 



figure 11). This apparent decrease in 7 might be due to a loss of numerical accuracy, but 
we are led to suspect the possibility that this is a true physical effect. This is suggested by 
the behavior of A. Let us analyze these two possibilities in more detail. 

What could be the numerical cause for an underestimation of 7? In our scheme the 



9 We do not address the details of this fitting since it has indicative aims only. Yet, one must bear in 
mind that the robustness of fitting is judged by the allowed variation of the fitting coefficients that still 
confines the fit within the errorbars. The proposed (quadratic) fit is significantly better in that sense than 
all other polynomial fits, making other behaviors than a minimum in A possible but less favored. 
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log(p.-p) 

Figure 8: The logarithm of the maximal curvature on the axis as a function of log(p* — p) in 
Because of small number of "wiggles" and their large amplitude, the linear fitting approach fails to 
accurately estimate the slope. In this case we employ the fact that the period of "wiggles" (~ 2.85) 
is equal to A/ (2^) and use this to find 7 ~ 0.44 from the more accurate A. This approach is seen 
in lower dimensions to be consistent with the regular way to obtain 7. 




-20 -18 -16 -14 -12 

log(p.-p) 



Figure 9: D = 8 Normalized plot of \og(R ma x) + 2j log(p* — p) (stars) and a sine wave (solid line). 
The period of "wiggles" is about 6.18 and it agrees well with theoretical prediction A/7 ~ 6.22. 

threshold amplitude p* can be obtained with somewhat decreasing accuracy as the dimen- 
sion grows. Unfortunately our evaluation of 7 relies heavily on the value of p*; the abscissa 
in figures ||, [?] and || is — log(p* —p)- Therefore, the quality of the linear fits in these figures 
and hence 7, depends acutely on how well we estimate p*. The effect is such that if the 
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Figure 10: A, including the errorbars, vs the spacetime dimension D. The solid line designates a 
2 nd degree polynomial fitting, the dotted lines show the 95% prediction bounds. The speculative 
extrapolation to higher dimensions suggests that A reaches a minimum at D ~ 12. 




Figure 11: 7 and 7 maS s, with the errorbars, as a function of the space-time dimension D. The 
dominant dimensional dependence of 7 maS s — (D — 3)7 is essentially linear increase due to the 
factor D — 3, which masks the rather subtle variation of 7. The solid line designates a 4th order 
polynomial fit based on the data points in the range 4 < D < 10 (the 7 value for D = 11 is not 
used). The dotted lines show the 95% prediction bounds. The extrapolation for larger D's indicates 
the existence of a maximum somewhere in between 10 < D < 12. 

accuracy of p* is lowered from 1 part in 10 14 to 1 part in 10 10 a few percents decrease in 
7 is induced. This might well be the reason we observe a downturn in j(D) at D = 11. 
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Having identified the possible cause for a numerical error in 7 at D = 11, we can try to 
estimate this value independently of the numerical calculation at this dimension. To this 
end, we fit a polynomial (this time of 4th order) to 7 using only data from 4 < D < 10. 
By extrapolating this fit beyond the last (D = 10) data point we again observe a change in 
the trend: 7 seems to reach a maximum at about D = 11 and then decreases. Remarkably, 
the numerically computed 7 value at D = 11 complies with this extrapolated behavior. We 
tend to interpret this as a hint that the observed downturn in ^y(D) reflects the physical 
behavior of the system. This, as well as the behavior of A(D), marks the appearance of 
interesting and non-trivial dimensional dependence associated with critical collapse. (One 
might call the dimension where this occurs a "critical dimension".) If so, this intriguing 
phenomena clearly deserves better understanding. 

While we stress that the allusion to the existence of a special dimension where the 
behavior of the the echoing/scaling exponents changes is based on extrapolation, the phe- 
nomenon itself is not entirely unexpected. The Choptuik scaling, albeit in somewhat 
different circumstances, is conjectured to appear [14] in the black-string-black-hole system 
which does exhibit critical dimensions 10 . In that system, one considers black objects in a 
higher-dimensional spacetime with one compact dimension, R- 0-2,1 x S 1 . The known solu- 
tions are divided to black string solutions (whose horizon wraps the compact direction and 
so has the S 1 x S D ~ 3 topology), and black hole solutions (with spherical, S D ~ 2 , horizon 
topology) localized on the circle. Consider the space of static solutions in this system. The 
different phases of solutions (black hole and black strings) are distinguished by an order 
parameter in this space, and the transition between these phases denotes a change in the 
topology of the solution. The local analysis [14] in the vicinity of the topology-changing 
(or merger) point shows that the geometry of the spacetime there should resemble a time- 
symmetric version of Choptuik's critical solution. Formally, the dilaton field gu plays the 
role of the scalar field in Choptuik's case, and the equations are the same. In both cases 
the problem is essentially two dimensional. 

The similarity between these two systems lends insight into what happens near the 
merger point: the geometry here, as speculated in [14|, is DSS and it resembles the Choptuik 
critical solution 11 . However, we would much prefer to have information on the dynamic 
solutions describing this phase transition. An attempt to trace the fate of an unstable 
string 12 in 5D was undertaken in [p3| , p4| , where it was shown that before the simulation 
crashes the black string becomes extremely non-uniform with a very long and thin neck. 
The numerical solution in [23, 24 1 did not find any fine structure of the critical behavior 
near the pinch-off point. We believe that an improved code will discover that the behavior 
close to the pinching is similar to what happens in the axisymmetric near-critical collapse. 
In the latter case, evidence was given in p5[ that a non-spherical mode appears causing 



10 The critical dimensions in this system are: (i) D mGr gcr = 10, above which the local geometry near the 
merger point argued to be cone-like, and below which this cone-like behavior is spontaneously broken^, 
and (ii) D(2d order) = 13 above which the phase transition becomes of second order J(| Q . 

n We note, however, that it is still not quite clear how to relate the Choptuik solution emerging in the 
collapse situation to its time-symmetric version studied in 

12 It is known that a uniform black string becomes unstable j5^] if "too thin" (Actually the relevant 
parameter is the ratio between the Schwarzschild radius of the string and the radius of the compact circle). 
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bifurcations of the axisymmetric clumps of collapsing matter, which is reminiscent of the 
situation expected in the black-string pinching. 

In summary: we have obtained A(Z?) and 7(D) for 4 < D < 11 and found clues to the 
existence of a critical dimension where the behavior of these functions qualitatively changes. 
However, clues are still only that. It is important and interesting to improve the numerics 
and discover what really happens beyond D = 11. If extrema are present in the functions 
A(-D) and "f{D), what are the trends of these functions for even higher dimensions? Will 
they be discontinuous at some dimension, indicating a phase transition? Alternatively, if 
these functions are continuous will they tend towards some constant values? Perhaps A 
blows up at a certain dimension signaling that the solution ceases to be DSS and becomes 
continuously self-similar (cone-like) beyond that dimension? We leave these questions open 
for future work, that will shed more light on this stimulating phenomenon. 
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